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ABSTRACT 

Context. The observed physical properties of outflows from low-mass sources put constraints on possible ejection mechanisms. 
Historically, these quantities have been derived from CO observations using ground-based observations. It is thus important to inves- 
tigate whether parameters such as momentum rate (thrust) and mechanical luminosity (power) are the same when different molecular 
tracers are used. 

Aims. We aim at determining the outflow momentum, dynamical time-scale, thrust, energy and power using CO and H2O as tracers 
of outflow activity. 

Methods. Within the framework of the WISH key program, three molecular outflows from Class sources have been mapped using 
the HlFl instrument aboard Herschel. We use these observations together with previously published H2 data to infer the physical 
properties of the outflows. We compare the physical properties derived here with previous estimates based on CO observations. 
Results. Inspection of the spatial distribution of H2O and H2 confirms that these molecules are co-spatial. The most prominent 
emission peaks in H2 coincide with strong H2O emission peaks and the estimated widths of the flows when using the two tracers are 
comparable. 

Conclusions. For the momentum rate and the mechanical luminosity, inferred values are independent of which tracer that is used, i.e., 
the values agree to within a factor of 4 and 3 respectively. 

Key words. ISM: individual objects: VLA 1623, L 1448, L 1157 - ISM: molecules - ISM: jets and outflows - stars: formation - 
stars: winds, outflows 
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C^l 1- Introduction Some of the physical properties that can be determined from 

-y-! ' these observations are of fundamental importance in the under- 

y ; The emission from CO is a widely used tracer of outflow activity standing of the star formation process. For example, the ob- 

^ The lowest rotational transitions emit photons in the millimetre served ratio between the momentum inputs of outflows and the 

^ ■ and sub-millimetre part of the electromagnetic spectrum mak- luminosities of the central sources sets constraints on the pos- 

; ing CO relatively straightforward to observe using ground-based sible ejection mechanisms (lLadalll985h . Using CO as a tracer, 

, telescopes. Combined with the fact that CO is abundant with various parameters of interest have therefore been deduced for 

■ respect to molecular hydrogen, this has led to extensive map- several molecular oiitflows in the past (see e.g.FAndre et al. 199^ 

Cl ' Pi"g campaigns, covering entire outflows, during the last two iBachiUer et al.ll990ll200lh . It should be noted, however, that un- 

C3 . decades. Rotational transitions of H2O and H2 are, on the con- certainties in the outflow mass and line-of-sight inclination angle 

trary, not easily observed from the ground and recent observa- can introduce large errors when estimating the total energy and 

tions have primarily relied on the use of space-based telescopes momentum of the outflowing material. 

(e.g. ISO, SWAS and Odin). However, the liinited achievable ^.^j^.^^ ^j^^ framework of the Water In Star-f onninf> re- 

spatial resolution on these missions hampered the possibility to .^^^ ^.^/^ ^^^^^j^^^ ^^^^^^ Ivan Dishoeck et all l20Tlb key 

interpret the data. The situat ion has improve d since the launch o f t^ree nearby molecular ou tflows from Cl ass 

Spuzer ^Werner etaI. |g004D and Herschel ^Pilbratt et al.| m). sources were ob served using the HIFI ^ de Graauw et aL l bOl l 

It IS now possible to observe shocked gas, whejetli^water abun- iR^elfsema et al.1 l20T2h instrument aboard Hers chel. The map 

dance is expected to be enhanced (see e.g. |Bergm et al. 19981 nh«Prv«tinn« nnv^r th. nntfln«,« T 11 ^7 (iTWmntnPt.l 



dance is expected to be e ntianced (see e.g. IBergiiLet at t W«t . observat ions cover the outflo ws LI 157 dUmemoto et al. 
I Kaufman &NeufeM | l996|) , at a much higher spatial resolution. f|^^ ^1448 (iBachiller et al.ll990l) and VLA 1623 (I Andre et al 



1199(1 119931) . Observations of the ground-state transition of 
* //er.sc7!e/ is an ESA space observatory with science instruments orf/zo-water, H2O (1 10 - loi), at 557 GHz have bee i i dis- 
provided by European-led Principallnvestigator consortia and with im- cussed, in conjunction with PACS dPogUtsch et alj |2010|) 
portant participation from NASA. observations of the H20(2i2-loi) line at 1670 GHz in 
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Table 1. Observations carried out with HIFI. 



3. Results 



Target 


O'lOOO 


(52000 




Obs. ID 




(hr:min:s) 


(deg:min:s) 


(sec) 




L 1448° 


03:25:38.40 


+30:44:06.0 


6468 


1342203200 


VLA 1623* 


16:26:26.38 


-24:24:31.0 


6552 


1342204010 


LI 157 


20:39:06.20 


+68:02:16.0 


6440 


1342210068 



Notes. (^^ ^Nisini et ai] ( 1201 3h iBierkeli et al] ( 1201 2h 



a series of paper s dNisini et all 20131 l iataiia et ai.[ izui j| 
iBjerkeli et all (20 121: ISantangelo et al.l 120121: IVasta et al.ll2012l) 



Tafalla et al 



2013t 



Th ese outflows were als o observed wi t h Spitze r and discussed 
in iNeufeld et al.' ('20091) iNisini et al.l (1201 Obi) . iGiannini et all 
(12011) and Bierkeli etatl ( |2012|) . The analysis of the H2 and 
H2O emission suggests that H2O correlates to a larger extent 
with H2 than with low-/ CO, both in terms of spatial distribu- 
tion and excitation conditions. The latter molecule is primarily a 
tracer of the cold, entrained ambient gas while H2 and H2O are 



probing the shocked gas at elevated temperature s (iTafalla et al 



120 131: INisini et al.ll2010al: ISantangelo et"ani2012l) . iBierkeU et al 
(I2OI2I) compare the physical parameters derived from CO and 



H2O for the VLA 1623 outflow. The authors conclude that the 
force driving the outflow seems to be independent of whether 
CO or H2O and H2 are used as tracers for the outflow. The analy- 
sis, however, was based on the observation of a single object and 
only the north-western outflow lobe of that object. In this letter, 
we therefore extend this analysis to include also the other out- 
flows that have been mapped with HIFI, i.e. L 1 157 and L 1448. 

We determine the outflow dynamical age, momentum, mo- 
mentum rate, energy and mechanical luminosity for LI 157, 
L1448 and VLA 1623. The physical properties derived from the 
observations of CO(l-O) and CO(2-l) are taken from previously 
published work. Using the combined spatial and kinematical in- 
formation obtained from the observations of H2O and H2 we 
then do the same analysis and compare the results. 



2. Observations 

The Spitzer observations of the purely rotational H? transitions, 
H2 - S(0) to S(7) were originally presented in iNeufeld et al.l 
(l2009h and we refer the interested reader to that publication. 

The mapping observations of the H2O (lio - loi) emission, 
from the outflows tar geted during the WISH program, have al- 
ready been published dNisini et al.ll2013t IBierkeli et al.ll2012h in 
recent papers, with one exception, viz. LI 157. All three maps 
were obtained with the HIFI instrument using the "On-The- 
Fly (OTF) Maps with Position-Switch Reference" observing 
mode. The maps cover the 5' by 2' regions closest to the cen- 
tral source of the outflows, which means that a large proportion 
of the L 1448 and L 1157 outflows are covered. In the case of 
VLA 1623, only a small part of the southeastern flow is covered 
by the observations and consequently we exclude this outflow 
lobe in the analysis. The data were obtained from the Herschel 
Science Archive (HSA, v8.2.1 of the pipeline) and standard data 
reduction methods such as baseline removal and production of 
maps were done in XS and Matlab. The unce rtainty attributed 
to the calibration is estimated to be ~10% dRoelfsema et al 



2012) and the main beam efficiency at 557 GHz is 0.76 dQlberg 
2OIOI) . The Herschel observations discussed in the present letter 
are summarised in Table [1] where the corresponding observation 
identification numbers are listed. 



For the work presented here, we rely on an estimate of the out- 
flow width and length (when estimating the mass and the dy- 
namical timescale) as measured from the H2O emission, which 
is why the HPBW of 38" at 557 GHz turns out to be a prob- 
lem. We therefore use a Statistical Image Deconvolution (SID) 
technique to im prove the spatial resolution of the HIFI data 
(lRvdbeckl2008b . The deconvolved HIFI maps are presented and 
overlaid on the Spitzer H2 emission in Fig. [1] From the decon- 
volution of the maps, we conclude that the widths of the flows 
are overestimated when using the original datasets. However, the 
apparent width of the northern flow in L 1 157 does not change 
by a significant amount when the SID method is used, and this 
part of the flow is likely broader than the southern part. Apart 
from this single outflow lobe, the outflow lobes have widths of 
the order -20". 

The spatial distributions of the H2 emission regions are in 
good agreement with the regions responsible for the H2O emis- 
sion. The most prominent peaks in the H2 emission, from all 
three outflows, are coincident with strong H2O peaks. However, 
we note that the H2O peaks towards the central sources of L 1448 
and VLA 1623 are not accompanied with strong H2 peaks. In 
the original H2O data, we detect a weak red-shifted component 
in the blue-shifted part of the L 1 157 outflow. After deconvolu- 
tion, this emission stands out and the brightness is significantly 
enhanced. It will, however, not be discussed further here. Note 
that the H2O line pr ofiles presented in Nisini et al. (2013) and 
IBierkeli et al.l d2012l) trace velocity ranges equivalent to those of 
CO. Hence, H2O does not seem to be special insofar as it traces 
gas moving at higher or lower speeds than CO. The maximum 
detected velocity for both molecules are in close agreement (i.e. 
within a few km s"'). 

4. Discussion 

In this section we compare the physical properties derived from 
CO observations to those derived from the H2O and H2 obser- 
vations. Although a resemblance (both what regards the spatial 
distribution and the excitation conditions) between H2O and H2 
is well established at this point, it should be noted that this anal- 
ysis relies on the assumption that also the kinematics are the 
same. Th e physical propertie s of the CO o utflows have been pre- 
sented in Andre et al. (1990; VLA 1 623), iBachiUer et all dl990l 
L1448) and Bachiller et alj d200lL LI 157). The results pre- 
sented in these papers are, however, corrected here, using the 
most recent method for determining these quantities. The ques- 
tion of how outflows are accelerated is still a matter of debate. It 
has been argued, however, that the highly collimated outflows 
that are seen to emanate from the very youngest objects (i.e. 
Clas s sources) are best exp l ained by jet-driv en models (see 
e.g. iMasson & CheminI 119931: lArce et al.ll2007l) . We therefore 
calculate the dynamical time-scale, momentum and kinetic en- 
er gy (t^, P and E resp ectively) in the same manner as suggested 
bv iDownes & CabritI d2007l based on numerical simulations of 
jet-driven outflows). The inclination angles of the flows with re- 
spect to the plane of the sky are assumed to be 21°, 15° and 9° 
for L 1448, VLA 1623 and L 1 157, respectiv ely dGirart & AcordI 
I2OOII: iDavis et alJll999l: iGueth et al.lll996h . It should be noted 
that uncertainties in the inferred inclination angles does not af- 
fect our main conclusions since the introduced errors are equiv- 
alent for CO and H2O. The largest relevant contribution to the 
error budget is instead attributed to the est imated molecular 
mass M of the outflows. Also, as noted by iDownes & CabritI 
d2007h . to use the intensity-weighted velocity <v>iobe (see e.g. 
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Fig. 1. Colormaps show the sum of the H2 S(0) - S(7) emission (normahzed to 1 x lO"-' erg cm"- s"' sr"') and contours of the 
H2O (1 10 - loi) emission are overlaid. Red contours are from 0.06 to 19 K km s \ 0.09 to 55 K km s ' and 0.05 to 7.9 K km s ' for 
L 1 157, L 1448 and VLA 1623, respectively. Blue contours for the same sources are from 0.06 to 69 K kms"', 0.09 to 23 K km s"', 
and 0.06 to 10 K km s ' . Offsets in the maps are with respect to the central sources, with coordinates tabulated in Table[T] The black 
dots indicate the readout positions for the OTF maps. 



lLada&Fichlll996h . averaged over the outflow lobe when deter- 
mining the dynamical time-scale may overestimate the age of 
the flow (unless the flow is inclined nearly 90° with respect to 
the line of sight). Therefore, and as discussed by the same au- 
thors, we instead estimate the outflow dynamical ages from the 
maximum velocity, i.e., dividing the extent of the flow with the 
maximum inclination corrected velocity - Liobe/>Jma.\- The 
maximum ve locity traced by CO is taken from the literature 
dAndre et al.i ri990: Bachill eret all [1990. 2001). When estimat- 
ing the momentum of the flow we use the intensity-weighted ve- 
locity and we do not apply any inclination correction. For the ki- 
netic energy we use a correction fa ctor of 5, i.e., the uncorr ected 
value is multiplied by this number (iDownes & Cabritll2007[ their 
Fig. 4 & 5). 

The originally published parameters based on CO data are 
corrected for most r ecent distance estimates, 250 pc to L 1 157 
dLoonev et al.ll2007h . 232 pc to L 1448 (Hirot a et al.ll201 ll) . and 
120 pc to VLA 1623 (iLombardi et al.il2008l) . The length of the 
flows are taken to be equal to the maximum extension of the 
half-power integrated intensity of each lobe. As mentioned ear- 
lier, the angular width of all lobes (except the northern flow of 
L 1 157) are close to 20". The molecular mass of the H2O emit- 
ting regions is calcu l ated from the e stimated H2 column densit y 
(lNisini etal.ll2010bt iGiannini etaT]l201ll: iBierkeU et aLll2012h . 
We assume the typical velocity of the wind to be 200 kms 
when estimating the mass-loss rate. All inferred quantities are 
summarised in Table |2] 

The highest detected velocity of the H2O emission line pro- 
files is consistent with previous estimates from CO observa- 
tions (see Tab. |2]i. The dynamical ages of the outflows are of 
the order 10^ years and this confirms that these sources are 
very young. They may therefore still be in the formation stage 
where material is accreted onto the proto-stellar condensation. 
The dynamical age can be compared to the shock propaga- 
tion time-scale (the time the gas is in a sho cked state) for var- 
ious types of plane-parallel shock models (iBergin et al.lll998t 
flower & Pineau des Foret^ l2010l) . For C-type shocks, time- 
scales are of the order 10^ years for a pre-shock density of 
10"* cm"^ and even shorter than this for J-type shocks or when 



13 - 



11 - 
10- 
9- 




M <^>,„fc ■»„„,, P E dP/dl dE/dt dM/dt 



Fig. 2. Ratio between parameters inferred from CO and H2O ob- 
servations. The markers indicate the ratios for L 1448 (squares), 
VLA 1623 (diamonds) and L 1 157 (circles). Red and blue colors 
are for the red and blue outflow lobes, respectively. Note that a 
few of the symbols are superpositioned. Dashed lines show the 
mean and boxes indicate the standard deviation. 



pre-shock densities are higher The dynamical ages can there- 
fore be longer than the time-scales for shock propagation. For 
the outflows discussed here, the water emitting regions trace gas 
at densities of ~10^ cm"^ and temperatures higher than 100 K . 
Under these conditions, the time-scale for freeze-o ut onto grains 
should be longer than the dy namical time-scale (IBergin et al.l 
119981: iHollenbach et al.ll2009l) . Thus, it is possible that once the 
H2O abundance is enhanced due to the presence of a shock it 
will not be significantly reduced in the post-shock region. This 
is in agreement with the fact that the H2O emission does not 
trace gas at particularly high velocities. The thrust dP/dt (mo- 
mentum rate) and power dE/dt (mechanical luminosity) derived 
from H2O and H2 are in close agreement with previous estimates 
based on CO observations. However, although the deviations are 
small, the thrust is on average 3 times higher and the power is 2 
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Table 2. Physical parameters derived from CO and H2O observations. 



Source: 

Molecular tracer: 
Red-shifted lobe 
Total mass (Mq): 
^lobe (pc) 
<i/>i„be (kms"'): 

(kms"'): 
fdyn (yr): 

Momentum {Mq km s"'): 
Kinetic energy (erg): 
Momentum rate (Mg km s^' yr" 
Mechanical luminosity (Lg): 
Wind mass-loss rate (Mg yr~'): 

Blue-shifted lobe 
Total mass (Mg): 

L\obe (pc) 

<t/>i„be (kms"'): 
i/ma, (kms"'): 
«dyn (yr): 

Momentum (Mg km s"'): 
Kinetic energy (erg): 
Momentum rate (Mg km s"' yr" 
Mechanical luminosity (Lg): 
Wind mass-loss rate (Mg yr"'): 



LI 448 



VLA 1623 



L1157 



CO 


H2O 


CO 


H2O 


CO 


H2O 


1.6 X 10"- 


6.3 X 10"-^ 


1.1 X 10"- 


1.6 X 10"^ 


8.4 X 10"- 


1.3 X 10"^ 


0.17 


0.17 


0.09 


0.09 


0.21 


0.22 


75 


64 


27 


60 


75 


152 


195 


218 


116 


116 


143 


192 


8.7 X 10- 


7.7 X 10- 


7.2 X 10- 


7.6 X 10- 


1.4 X 10-3 


1.1 X 10' 


0.42 


0.14 


0.08 


0.03 


0.98 


0.30 


5.6 X 10''^ 


1.6 X lO"" 


2.7 X 10"-^ 


2.0 X 10"-'' 


5.7 X 10"" 


3.5 x 10"" 


4.8 X lO""* 


1.9 X 10"" 


1.1 X 10"" 


3.3 X 10"^ 


7.0 X 10"" 


2.7 X 10"" 


5.3 X 10" 


1.8 X 10" 


3.1 X 10"' 


2.1 X 10"' 


3.4 X 10" 


2.6 x 10° 


2.4 X lO"** 


9.4 X 10"' 


5.5 X 10"' 


1.7 X 10"' 


3.5 X 10"" 


1.3 X 10"" 


2.0 X 10"- 


2.8 X 10-3 






1.2 X 10-' 


9.4 x 10"-3 


0.16 


0.08 






0.12 


0.12 


-75 


-94 






-29 


-49 


-195 


-195 






-68 


-128 


7.8 X 10^ 


3.8x10- 






1.8 X 10-3 


9.4 X 10^ 


0.53 


0.09 






0.54 


0.07 


7.1 X lO^-* 


1.6 X 10"" 






1.2 X 10"" 


2.7x 10"3 


6.7 X 10"" 


2.5 X 10"" 






3.0 X 10"" 


7.6 x 10"' 


7.5 X 10" 


3.4 X 10" 






5.8 X 10"' 


2.4 x 10"' 


3.4 X 10"" 


1.2 X 10"" 






1.5x10-" 


3.8 x 10"' 



Notes. Inclination angles with respect to the plane of the sky are 21°, 15°, and 9° for L 1448, VLA 1623 and L 1 157, respectively. The velocity of 
the wind is assumed to be 200 km s"'. 



times higher when deduced from CO observations (see Fig. |2). 
The reason for this is not entirely clear but it may be that the 
mass of the H2O emitting gas in reality is slightly larger than 
the values estimated from this work. The wind mass-loss rates 
are estimated to be -10"'' Moyr"' for L 1448 and L 1157, and 
one order of magnitude lower for VLA 1623. Assuming that the 
luminosity of the central sourceQ is entirely due to accretion and 
adopting a radius to mass ratio of 5, this leads to simila r values 
for the mass-accretion rates (see e.g. IStahler et al.lll980l) . 

5. Conclusions 

The inferred values for power agrees to within a factor of 3, 
when using the different molecular tracers. Similarly, the esti- 
mated values for thrust are equal to within a factor of 4. This 
may indicate that the ejection mechanism that is responsible for 
the motion of the cold gas, is the same as the one that sets the 
warm gas into motion. It is, however, also possible that this is 
a characteristic of young outflows. This study also reveals that 
the emission from H2O traces a gas component that presumably 
operates on the same time-scales as CO, i.e. H2O is not detected 
at very different velocities compared to CO. The fact that the 
estimated values are the same (within a factor of a few) for all 
observed sources supports previous estimates based on CO ob- 
servations. We also find it likely, that ground-based CO obser- 
vations are adequate, when assessing the impact of outflows on 
their environment. 
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